Growth and Decay of Magnetic Fields 
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ABSTRACT 


We investigate the magnetohydrodynamic properties 
of random, stationary turbulence. If the fluid is highly 
conducting and the magnetic field weak, the logarithm of 
the field strength B is shown to be a random-walk variable 
to which the central limit theorem applies. The expecta- 
tion value of B is then shown to increase exponentially 
with time. The behavior of infinitesimal material dis- 
placements which move frozen in the fluid, and which are, 
on the average, stretched by ' the turbulence is investi- 
gated. It is well-known that for high conductivity the 
quantity J3/density in a given fluid element remains pro- 
portional to the length of such a displacement chosen 
initially parallel to B. Further, we show, with the use 
of a special gauge condition, that the magnitude of the 
vector potential A/density remains proportional to the 
area of an infinitesimal material parallelogram chosen 
perpendicular to A. The central limit theorem for two 
random variables then gives the joint probability density 
for A and B. If B has no large-scale Fourier component 
(imposed by boundary conditions) , the scale of the field 
is A/B. The resulting ohmic dissipation grows much 
faster than the random-walk energy input from the fluid. 
Plausible extension to high ohmic dissipation is made, and 
we conclude that in the absence of a large-scale Fourier 
component, the ohmic loss destroys the field. An experi- 
ment is proposed to test these conclusions. 
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I . INTRODUCTION 


The problem of the amplification and maintenance of 
magnetic fields in turbulent conducting fluids has long 
interested hydrodynamicists and astrophysicists. One 
opinion about the behavior of such systems that is gen- 
erally agreed upon is that under suitable conditions an 
initial magnetic "seed" field will, on the average, be 

1 

amplified by the turbulent motion (Biermann and Schluter , 

2 3 *+ 5 6 

Batchelor ' , Syrovatski , Saffman ' , Kraichnan and 

7 8 

Nagarajan , and Pao ) . However, the amplification process 
is also accompanied by a steady decrease in the scale of 
the field, so that the ohmic losses may in the long run 
overtake the amplification and wash the magnetic field out 
of existence. 

Many authors have conjectured that turbulence will 

indeed finally destroy a magnetic field (Zeldovich 9 , 

Saffman 5 , and Parker 10 ), but no convincing proof of this 

9 

conjecture has appeared. Zeldovitch proved that two- 
dimensional turbulence will, under certain rather restric- 
tive boundary conditions, destroy a magnetic field, but it 
seems impossible to generalize his method to three dimen- 
sions . 

The present paper is concerned with setting up a 
statistical formalism for stationary random turbulence which 
is exact in the case of high conductivity and weak magnetic 


fields . 



By "weak" we mean that the Lorentz forces may be 

neglected in the fluid equations of motion. We use the 

central limit theorem to derive an exponential law of 

magnetic field increase with time by showing that the 

logarithm of the field strength in a given fluid element 

is a good random-walk variable. Biermann and Schluter 

were the first to surmise the exponential increase of the 

field strength. Our analysis is related in some ways to 
10 

that of Parker , who also derived an exponential increase, 
and we are able to provide a more nearly rigorous basis for 
some of his results, with which we are in substantial agree 
ment. 

We derive a similar exponential law for the magnitude 
of the vector potential and are then able to show that the 
ohmic dissipation increases at a much faster rate than the 
magnetic field energy. Whether or not a final steady state 
can be maintained is shown to depend in an important way on 
boundary conditions, and we argue plausibly, although not 
rigorously, that if the magnetic field is not imposed from 
the outside by boundary conditions, the turbulence must 
eventually destroy the field, no matter how large the con- 
ductivity. 

If the field is initially strong enough to suppress 
the turbulence to some extent, this conclusion may not 
necessarily be valid. 


- 3 - 



A feasible laboratory experiment using liquid sodium 
is proposed which might show whether or not the general 
conclusions of this paper are valid. It would also pro- 
vide a test of Batchelor's criterion for the maintenance 
of turbulent fields, with which the present paper disagrees. 

2. LOGARITHMIC FIELD STRENGTH AND VECTOR POTENTIAL 
AS RANDOM WALK VARIABLES 


The equations of motion of the magnetic field B^ are, 
in Gaussian electromagnetic units, 


dB _ 

dt 


+ (U.V)B = (B.V)u - B(V.U) - VX(XVXB) , 


( 1 ) 


where u is the fluid velocity field and X — (4na) is 
the magnetic diffusivity, a being the conductivity in 
sec/cm 2 . If we assume that the conductivity is infinite, 
then the equations become, in three-dimensional Cartesian 
tensor form, with a comma indicating partial differentiation. 


dB 

dt~~ = B n U a,n “ B a u n,n ” T an B n' 


where T flb h u a , b - S ab u n , n - 


(2) 
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The material derivative d/dt follows the motion of a 
given fluid element as it experiences the various velocity 
shears and compressions. In what follows we investigate 
the behavior of the field and vector potential in the 
given fluid element as it moves about in the turbulent 
medium. Thus we suppress the spatial dependence of the 
dependent variables and consider explicitly only their 
time dependence. 


2.1 The Field Strength 

We consider the fluid motions to be turbulent and 
random, and thus we assume that the tensor T^(t) is a 
random function and that it is, for weak fields, inde- 
pendent of the field strength B(t) in the fluid element. 

Let us substitute = Bn in Eqn. (2) , where n (t) is a 

a a a. 

unit vector, and then form the scalar product with n a , 

using the fact that n dn /dt = 0, to obtain 
J cl a . 


d In B 
dt 


n a n b T ab' 


(3) 


It is thus apparent that lnB(t) can be used as a 
random-walk variable, since the right-hand side of this 
equation depends only on and on the direction of 

the field. We now break up the time axis into steps of 
constant length fit, such that the change in InB over fit 
is essentially uncorrelated with the changes over the 
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previous steps. If t is the time scale of the larger com- 
ponents of the turbulence, we might take St 2 t , It is 
clear that since T at) (t) should be continuous in t, we can 
never completely lose the correlation with previous values 

of T a b* However, for St >, 2x, these correlations should 

1 1 

become very small 

We can now use the central limit theorem to obtain 

the probability density for InB s b after a large number 

N of time steps. Let S^b = S^lnB be the change in InB 

over the jth time interval. For stationary turbulence, 

the probability density function for S^b will not depend 

on j, and we thus write y e ^Vb') and 0 ^ 2 e ^(Sjb) 2 >- y 2 >0, 

where y and 0 1 do not depend on j. Under these circumstances, 

1 2 

the central limit theorem states that the probability 
density for the sum Ab = Ij =1 6jb = b (N) - b(0) for large N 
converges to the normal density 

f(Ab;N) = (2Tr9 1 2 N)“ Js exp [- (Ab-Ny ) 2 /20 : 2 N J . 

We now wish to argue that y = 0 . Since the turbulence 
is random, the components of the tensor should have zero 

expectation value, and since in any case all "memory" of the 
initial field direction n a ° is very quickly lost, it seems 
impossible even in the case of anisotropic turbulence that 
the right-hand side of Eqn. (3) should show any statistical 
tendency to be either positive or negative. Thus we 
sen y = 0, and the normal density for Ab simplifies to 
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f(Ab;N) - (2ir9 1 2 N)’ 33 exp [ ~ ( Ab) 2 /20 1 2 N ] , 


(4) 


3 

This is at variance with the assumption of Batchelor, 
that y > 0 . 

We may now derive the exponential law of increase 
for the expectation value of the field strength. If B 

o 

is the initial field strength in a given fluid element, 
then after N time steps, the expectation value of B in 
the same fluid element is = B q and we easily 

obtain 



B ( 2 TT 0 
o 



+oo 

dx exp 

— 00 




= B o exp(^e 



(5) 


Since N = t/6t, we see that ^B^ increases exponen- 
tially with the time. 

What is the order of magnitude of 9j? If we take 
<5 1 = 2 t, then probably 0 1 £ 1, the exact value depending 
on the type of turbulence. 

It is easy to find the law of increase for the 
magnetic energy density. We write 


<TO- sf- <' b2 >= %<> 24b > c o e X p(2 9l 2 N). (6) 

Thus <(b 2 >> (<B» 2 . 

Note, by contrast, that the conclusion y = 0 
implies that the expectation value of the change of 
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InB vanishes: ^AlnB^ = = 0. Thus, in a given fluid 

element, the field strength is just as apt to decrease as 
to increase. However, given that an increase in B occurs, 
the expected change in |b| is much larger than when a 
decrease occurs. 


2.2 The Vector Potential 

Having found the normal distribution for InB, we now 
do the same for the logarithm of the vector potential mag- 
nitude A, again under the assumption of infinite conductiv- 
ity. We first show that with the choice of a special gauge 
condition the vector potential A in a given fluid element 
remains perpendicular to an infinitesimal parallelogram 
formed by two material line elements which move frozen in 
the fluid element. l.e., if we define the parallelogram 
by two infinitesimal material displacements Sy and <$£ such 

that initially A c< <5y / <$z, then this proportionality holds 

***** *~ 

for all time. Furthermore, the quantity A/P remains pro- 
portional to the area of the parallelogram. 

The special gauge condition chosen is written ja .A. = c\p , 
where c is the velocity of light and is the scalar poten- 
tial. Writing Ohm's law as a = E + c ^u X taking 
a -> co, and using E = - 3A/3ct - and the above gauge 
condition, one finds 

§& = - (A.V)U - A>((vXu) = - <vu) .A. 

Qt v* *** ******* 
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1 3 

However it is easily shown that the equation of 
motion of the quantity p<5S = p6y/$z is identical to the 
above equation for A. Therefore A/p remains proportional 
to the infinitesimal surface element <5S, provided that the 
above gauge condition is employed. 

We may immediately recover the conclusion of 
9 

Zeldovitch that for two-dimensional incompressible motion, 

A_ is not changed by the motion itself, but only by the 
ohmic dissipation. In his analysis the only important 
component of is the one perpendicular to the plane of 
the motion, and since the motion is supposed incompressible, 
the areas of the parallelograms remain constant, and by our 
analysis A then remains constant. 

We now derive the probability density for the area of 
the parallelogram in the three-dimensional case. Since the 
parallelogram is infinitesimal, the shears and strains become, 
in the limit 6S-*0, spatially uniform over its surface, so that 
it remains a parallelogram. In Figure 1 we have indicated how 
the parallelogram may change over the time interval <5t, and we 
define three parameters which specify this change: (1) A homo- 

geneous horizontal scale change, which transforms AB-*A ' B (2) 
a homogeneous vertical scale change, which transforms CE+C ' E ' ; 
and (3) a pure transverse shear which changes SB’+A'D', but 
does not change AB or CE. Of course, the orientation of the 
parallelogram changes as well, but this is of no interest here. 
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The area of the parallelogram is S = AB.cl - , 

S' = A r B' .C'^' , and thus the transverse shear does not 
change the area, and we may ignore it. We are left with 
the two perpendicular scale changes 6a and 60, which 
change InS additively over the time interval 6t; i.e., 
a given scale change (6a, 6 0) always induces the same 
6lnS, independently of the initial value of S. Then de- 
fining 6a and 60 such that each represents the change of 
the In of a unit length, we set 6lnS = 6a + 60, and there- 
fore 6 In (A/p) = 6a + 60. Now, for stationary turbulence, 
the density may fluctuate, but in the long run we expect 
p # constant, since 61np is negatively correlated with 
previous 61np. Thus we ignore the density fluctuations 
and treat the fluid as incompressible. Therefore, we con- 
clude finally, for our purposes, 

6lnA £ 6a + 60. (7) 

l k 

However, it is well-known that the quantity B/p 
is proportional to the length of an infinitesimal 
displacement parallel to which moves frozen in the 
fluid. Therefore, since we have set p~ constant, the 
6b introduced in Sec. 2.1 is the same type of param- 
eter as the 6a and 60 defined directly above, since it is 
likewise the change in the In of the length of a material 
displacement. 

Therefore, except for the fact that 6a and 60 are 
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mutually perpendicular, the parameters fia, 50, and fib are 
statistically identical. Hence ^(60) 2 ^ «= <^(sb ) 2 ^ 

O' 

= 0 . Also, (&a,y = ^50^ = ^fib^ « 0, as in Sec, 2.1, 

2 

We define the correlation P x = </<Sa6 0^/O 1 , where -1 < Pj < + 1. 

Following the procedure in Sec. 2.1, we can now use a 

1 5 

two-dimensional form of the central limit theorem to obtain 

the joint probability density of the sum of N >> 1 two- 

dimensional variables ( Aa , A 0 ) = (<SjC*,5j8). Since the 

joint single-step density for (6^a,6j0) does not depend on 

previous values of these variables, the joint density for 

(Ac*, A 3) converges for large N to the bivariate normal density, 

1 6 

which for this case simplifies to 


f ( Act , A 0 ;N) 


= 2ir6 1 2 N(l - P i 2 )* 5 “ 1 X 


exp 



2 N(l-p 1 2 )]“ 1 [(Aa) 2 -2p 



( 8 ) 


What is the sign of p 1 ? Since the volume of the fluid 
element is approximately conserved, we expect fi^a and 5^0 to 
have different signs statistically, and hence p is negative. 

If we assume that the correlation is caused only by 
volume conservation, we can derive an exact value for p^. 

Let fi § be the third perpendicular logarithmic scale change. 
Since pure shears do not change the volume of the fluid 
element, the joint probability density f (Ac* , A0 , A § ;N) 
contains a delta-function factor 5 (Aa+A0+A§) ; and since 
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the arguments must appear symmetrically , we can write 


f (Aa,A3,A§;N)=C exp ^-D^(Aa) 2 + (AS) 2 + (AS) 2 J jfi (Aa+A3+A§) . 


We then integrate over A§ to get the marginal density 


f(Aa,A3;N) * 


dA§f (Aa,A3,A§ ;N) 


58 C exp j-2D (Aa) 2 +AaA3+ (A3) 2 


! ]}- 


Comparison with Egn. (8) gives the result 


Pi = ~h • (9) 

At this point we can make contact with some of the 
1 0 

results of Parker , who also derived the exponential 
increase of the magnetic field by considering the stretch- 
ing of material line elements. He set up a Fokker-Planck 
equation by supposing that the relative length change of 
a line element over a time step (large enough for the 
correlation with adjacent time steps to be small) was very 
small, but he recognized this assumption to be erroneous. 

Our use of the central limit theorem does not involve such 
an assumption, and also has the advantage of treating the 
local perpendicular scale changes by means of a joint 
probability density, instead of having to discuss the ap- 
proximate evolution of a macroscopic flux rope. 
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It is useful to derive an expression for the prob- 
ability density of Aa = Aa + A3. One may change variables 
in the (Aa,A3) plane to Aa = Aa + A3, Aa' = Aa and then 
integrate over Aa», the Jacobian of the transformation 
being unity, to get 

f(Aa;N) = 4tt 0 i ^N (1+p i ) ""^expf- ( Aa) ^ 46i^N(l+pi) 

. J v L 1 ) 

2 2 

which is the normal density with variance 20 ^ N(1 +pi)= 02 N. 
If Eqn. (9) holds, then 0 2 = 6i« 

Since the magnitude A of the vector potential in a 
given fluid element is proportional to the area S of the 
perpendicular parallelogram, we have A/A q = S/S q , and 
therefore, from Eq. (7), AlnA = Aa + A3 = Aa. Hence 

<A>=A o <e ia >' 01 

A Q exp (Js6 2 Z N) , which is similar to 
Eqn. (5) for the field strength. 

We have used the same assumptions here as those 
employed in Sec. 2.1; namely, that both the dynamic re- 
action of the field on the fluid and the ohmic dissipation 
may be neglected. 

3. THE JOINT DISTRIBUTION FOR A AND B, 

AND THE OHMIC DISSIPATION 

In this section we bring together the results of the 
previous section in order to gain information about the 
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rate of ohmic dissipation of the field. The previous re- 
sults depend on the assumption that the ohmic losses are 
negligible , but we will be able to extrapolate, with some 
reservations, to the regime where the ohmic losses are 
considerable. The conclusion will be that, except when 
the boundary conditions force a permanent large-scale 
Fourier component of the field, the field must ultimately 
be destroyed. 

3.1 The Local Scale and the Correlations 
Between InA and InB 

We begin by finding an expression for the local 
scale i of the magnetic field. Dimensionally, the re- 
lation A £ IB suggests itself. Now, this l is really 
the scale of A, but if B possesses no large-scale Fourier 
component, so that the polarity of B^ alternates randomly 
in space, then l Z A/B is, on the average, the scale of 
B as well. Thus we may state that for each (r,t) , there 
exists a point r + Ar^r ,t) , such that 

A(r,t )B (r,t) = c o A(r + Ar,t) , 

where c is a constant geometric factor of order unity, 
o 

It may happen that at a particular (r,t) , the 
field is flat, so that l ■* » there. Then Ar will not 


V 


( 10 ) 


- 14 - 


exist. However, we interpret % as an average scale over a 
finite region, and since we are interested only in average 
quantities in what follows, our conclusions will still be 
valid. It is important to emphasize that if has a large- 
scale Fourier component, Eqn. (10) is not valid. In Sec. 3.3 
we present a counterexample to results derived from Eqn. (10) 
and show how other aspects of our reasoning may break down as well. 

Let us now discuss the correlation between A and B. Since 
A*B/p oc p <5x • <5y X 6z, which is the total conserved mass in the 
material volume element defined by the three infinitesimals, 
the quantity A'B. is time-independent except for fluctuations 
of the density, which as we have said we may ignore for sta- 
tionary turbulence. Thus A*B^ = AB cose ^constant in a given 
fluid element, where 0 is the angle between A and ]3. There- 
fore, except for the anomalous case _A-B ==• 0, it follows from 
InA + InB + In |cosd| = In | constant | that, for a given change 
of cos6 , InA and InB will be negatively correlated, and we may 
use the generalized central limit theorem as in Sec. 2.2 to 
write the joint probability density, after N >> 1 time steps, 
of a I nA = Aa and AlnB = Ab as the bivariate normal density 


f(Aa,AbjN) - . 2tt9 1 Q 2N (1-p 2^) ^ ^ X 


exp 


1 - 


2N (1-p 2 ) 


-1 


(Aa) 2 2p 2 AaAb (Ab) 

1 1 1 *"* J r " 1 ~ ' a 1 • " ' ' n 


2n 


L 02' 


0 2 0 i 


(ID 


0l‘ 


where 



26 i 2 (1+p \ ) 


and - 1 < p 2 < 0 . 


- 15 - 



3.2 The Ohmic Dissipation Rate 


Using the expression for l provided by Eqn. (10),. we 

now find the expectation value of the ohmic loss rate. At 

3 

a point £ the loss of energy in erg/cm sec is then 

,X_ t 2 . . ^ XB (r ) 2 _ XB(r ) 4 

A ^ O O O * 

41T "* 4tt Jl (r) z 4uc ^A (r+Ar ) ^ 

a* O 

Now in order to find the expectation value of this 
quantity, we must weaken the correlation between Aa and Ab 
somewhat, since they are taken a distance Ar_ apart. Thus 
in Eqn. (11), we substitute p 2 p 3 , where P 2 < P 3 < °/ 
and obtain, for a given fluid element, 



X 

4ir 




— 8p 3 Q i 9 2 + 8 0 2 


2 



( 12 ) 



But Eqn. ( 6 ) implies that the rate of energy density 

increase from the random walk is, on the average, d ^e^/dt 

= 2 e 8i 2 (6t)" 1 exp(29 1 2 N) . Since p 3 < 0, Eqn. (13) tells 
0 

us chat no matter how small J Q is, the ohmic losses even - 

tually become comparable to the random walk gain, and the 

field ceases to grow. In particular, if pi = ~H and Si =62 

(Eqn. (9)) , .then <J 2 >= J Q 2 exp 0 X 2 N (10-8p 3 )J , which grows 

2 

much faster than exp (20 1 N) . 
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May we now conclude from this that the field even- 
tually destroys itself? Not yet, for the probability 
densities derived for A and B depended on the assumption 
that the ohmic losses were negligible. Indeed, one pos- 
sibility is that the ohmic dissipation can serve to in- 
crease the seal© of the field, and thus an equilibrium 
state might be possible. This possibility is discussed 
more thoroughly in the next subsection in connection with 
a counter-example, which also entails the possibility that 
the reasoning leading to Eqn. (10) breaks down. 

It is reasonable to conclude, however, that if the 
ohmic losses do not counter the tendency for the scale to 
decrease, and if Eqn. (10) holds, then the field must be 
destroyed. The probabilistic analysis given above casts 
considerable light on the action of the turbulent fluid 
motion on the field, and it is fair to say that the random 
walk process which tends to increase the field strength 
must at the same time inexorably be accompanied by a more 
drastic increase in dissipation: If O ) increases, 

then B ^/ %" 'y must increase even faster, since the dis- 

-2 

persion associated with the A part of the probability 
density is increasing at the same time. 
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3.3 A Counterexample 


We now discuss a special case, in which a large-scale 
Fourier component is forced upon the field by boundary 
conditions, and we show how the preceding analysis breaks 
down. 

Consider the situation depicted in Fig. 2, where the 
fluid is broken up into three regions: regions I and III, 

where X = 0 , u = 0 , and region II, where X > 0 , u ^ 0 . Let 
E be a part of the boundary surface between I and II such 
that the magnetic flux ^ = f f ^ B^d£> 0. 

Then $ ^ is constant in time, and the flux conserva- 
tion law V ■ B = 0 implies that the flux <»> through the 
"cap" surface ^2 in region II is the same as Therefore, 

no matter what the nature of the turbulence in region II, 
and no matter how low the conductivity is there, the mag- 
netic field in region II must persist. Thus an equilibrium 

1 7 

state is possible. 

1 .8 . 

Machine integrations by Weiss suggest that m such 
a system the magnetic field may be expelled out of the 
main body of the fluid into a boundary layer. However, 
Weiss dealt with time-dependent velocity fields, and one 
can imagine that in the random turbulent case the expulsion 
tendency might be countered by turbulent transport of flux 
back out of the boundary layer into the main body of 
region II. 
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Where does the analysis in Sec, 3.2 break down in 
this case? As mentioned before, there is the possibility 
that the ohmic dissipation tends to keep the scale of the 
field from decreasing. This can be shown to be related to 
the fact that in this counterexample the field has a perma- 
nent large-scale component. 

Consider the very simple case of a diffusion equation 
2 2 

SB/3 t = 3 B/3x . If the field is purely sinusoidal, then 
B = B i e"" <t sin (x/ic) is a solution, and the scale of the 
field Z = B/(3B/3x) = (k) ^tan (x/iT) is time-independent. 

Now suppose the field has a large-scale constant Fourier 
component, so that B = B^+B^ sin (x/k) is a solution. 

Then the scale becomes i = B(t)/|B 1 e Kt cos (x/<) /<J which 
diverges as t-*«° . We see that, at least in this crude 
example, the diffusion term leads to an increase in scale 
if there is a large-scale component. 

There is also another possibility that we must reckon 
with, which may also operate to invalidate the conclusions 
of Sec. 3.2, and which deals with the assumption of Eqn. 

(10) and its application to the ohmic loss estimate. In 
writing Eqn. (10), we assumed, in essence, that the field 
has a random character such that the contribution to 
coming from the field in distant volume elements was very 
small. However, if B has a large-scale component super- 
imposed on a random, alternating-polarity component, then 
the analysis might break down because Ar^ could conceivably 
introduce unexpected positive correlations betw’een InA and InB. 
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It is not clear which of these two effects plays the 
greater role in the above counterexample. 

3.4 Extension to High Dissipation and Strong Fields 

Let us now go back to the situation in which there 
is no permanent large-scale Fourier component of the field, 
and for which the ohmic losses presumably do not operate 
to increase the scale of the field. 

As the ohmic losses grow, the field is no longer 
frozen in the fluid, and begins to "slip" relative to the 
fluid. Thus the random-walk build-up and the scale decrease 
both become less efficient. However, the general prob- 
abilistic analysis for the random-walk transfer of energy 

from the fluid to the field and for the ohmic dissipation 

2 

should still hold. Presumably 0i then decreases, and J q 
must be continually renormalized to larger values . But with 
each random-walk extraction of energy from the fluid, the 
scale decreases. Thus the energy is renewable up to a 
point, but in the absence of a large-scale Fourier compo- 
nent the scale of the field seems doomed to a continual 
decrease, until the field is extinguished. 

We may also try to extend the reasoning to the case 
where the field strength is high enough to control the 
motion of the fluid. Where there is no large-scale compo- 
nent of the field, the same conclusion may result: It may 
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be impossible to put random-walk energy into the field with- 
out a continual decrease in scale. However, since the mag- 
netic force term in the fluid equations of motion is 
% (4n Si) , the field is more amenable to distortion where 

l is large, and hence large fields and large scales may be- 
come positively correlated. It is not evident that this 
tendency can counter the general reduction of scale implicit 
in the random-walk process. 

4. A NECESSARY CONDITION FOR TURBULENT AMPLIFICATION, 
AND A POSSIBLE EXPERIMENT 


The conclusions of the preceding sections are, as we 
have said, by no means rigorously derived. Therefore, it 
would be very useful to devise an experiment to test 
whether or not they are valid. In order to do this, we 
first show a necessary condition, elucidated also by 
Syrovatski** , for initial amplification of a magnetic field 
by turbulence. It is then shown that a technologically 
feasible experiment using liquid sodium as the conducting 
medium may be of considerable interest in testing our 
results . 

Suppose that we are given a medium of conductivity a , 


in which there exists turbulence of root-mean-square 


velocity u and characteristic scale defined perhaps 


1 9 

by the expression 
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is the turbulent energy spectrum, and k is the wave number. 

is thus the scale associated with the energy-bearing part 
of the spectrum and is presumably the scale most significant 
in distorting the magnetic field. 

Now, any magnetic field carried with the fluid must, 
if it is to be amplified or maintained by the turbulence, 
have a scale which is less than or equal to the turbulent 
scale. Further, Eqn. (1) shows that for amplification we 

must have I (B^£) u| >I^X(XV)(B) | , or in approximate form 

2 2 
u/i. t y X/J l . But since a £ we must have u/* t > X/A t , or 

4Trou2, t > 1. (13) 

2 

Batchelor used the analogy between vorticity and the 
magnetic field to derive a necessary and sufficient condi- 
tion for turbulent amplification and maintenance of the 
field. His condition is that the kinematic viscosity v 
be greater than X, or 4irov > 1. The physical reasoning be- 
hind this inequality is as follows: The smaller is v, the 

smaller is the scale of the turbulence. The scale of the field 
is inherited from the scale of the turbulence, and thus if v is 
small enough the ohmic losses will overwhelm the growth rate. 
However, our analysis shows that if the field has no large- 
scale Fourier component, the scale of the field in any case 
becomes ultimately so small that the field is extinguished. 

Let us now see under what circumstances we may satisfy 
Eqn. (13) , which is at least a necessary condition for turr 
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bulent amplification. One of the most highly conducting 

fluids available in the laboratory is liquid sodium, for 

which 2 °a & 10 -4 sec/cm 2 . If we can generate turbulence 

for which = 1 cm, then Eqn. (13) will be satisfied 
3 

for u > 10 cm/sec. 

Suppose that the liquid sodium be contained in a 
rotating cylindrical vessel of radius 16 cm and a counter- 
rotating set of paddles be situated inside to stir up the 
liquid. The scale of the energy-containing eddies is 
presumably inherited (but somewhat reduced) from the 
boundary scales. Thus in this case 1 cm. If the 

container and the paddles are counter-rotated each at a 
speed of 100 rev/sec, then the boundary speed of each is 
u ~ 10 4 * * cm/sec. This should insure the satisfaction of 
Eqn. (13) , but higher speeds might of course be used. 

The magnetic amplification properties of the ar- 
rangement might be tested with various externally im- 
pressed magnetic fields, but it would be difficult to 
set up the situation described in the counterexample, where 
the external conductivity must be much higher than in the 
turbulent fluid. 

The viscous dissipation in the fluid would amount to 

, o .. 21 

roughly vu' /*- t erg/gm sec. For liquid sodium , 

3 

v ~ 10"* 2 cm 2 /sec, and if we assume 10 gm of material, 

the total viscous power is then 10^ erg/sec = 100 watts. 

The turbulent kinetic energy density is 

ispu 2 £ 5 X 10 7 erg/cm 3 , and if equipartion with the fluid 
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is ever reached, the field strength would amount to 

. 4 

B 3 * 10 gauss. The associated ohmic dissipation is then 
2 2 *"1 3 

XB (4ir* t ) erg/cm sec, and the total ohmic power for 10 3 cm 3 

13 6 

is 10 erg/sec » 10 watts. Clearly, it would be impossible 
to reach this high power rate required for equipartion. If, 
however, turbulent maintenance of the field is possible, con- 
trary to the results of this paper, then an average mechanical 
power input of 200 watts should maintain a field of 300 gauss. 

One might also test the possibility that a field which is 
initially strong enough to control the turbulence could be 
maintained. This possibility is discussed at the end of the 
previous section. Since large-scale fields of 3 X 10 4 gauss 
are not easy to produce in the laboratory, a more tractable 
combination of u and might be tried. 

A negative result for this series of experiments would 
substantiate the conclusions of this paper, but would not show 
that Batchelor's criterion is wrong, since 4ircrv £ 10~ 5 . How- 
ever, a positive result would show that both Batchelor's cri- 
terion and our conclusions are incorrect, and that turbulent 
maintenance is possible. 
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FIG. 2. The counterexample to the results of Sec. 3. 
A large-scale permanent Fourier component is forced on the 
field in region II by boundary conditions. 



